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ABSTRACT 

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the 
Herschel/HIFI instrument. No strong lines are detected down to 3cr levels in 0.5 km s"' channels of 4.2 mK for the 1 lo-loi line and 12.6 mK for 
the lii-Ooo line. We report a very tentative detection, however, of the lio-loi line in the Wide Band Spectrometer, with a strength of r„i, = 2.7 mK, 
a width of 5.6 km s"' and an integrated intensity of 16.0 mK km s"' . The latter constitutes a 6cr detection. Regardless of the reality of this tentative 
detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated 
regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane. 
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1. Introduction 

Because of its association with biology on Earth, water is one 
of the most important molecules in the solar system and beyond. 
However, the origin of water on Earth is highly uncertain. What 
is clear is that the distribution of hydrated rocks in the solar sys- 
tem ( Abe et aT]|2000[ l suggests that water resided in the vapor 
phase in the warm (Tdust > 100 K) inner solar nebula and is pre- 
dominantly condensed in the form of ice beyond the so-called 
snow -line between 2-2.5 AU ( |Hayashi|1981[|Abe et al.pOOO] !. 

The observation of cool water vapor in protoplanetary disks 
is hampered by atmospheric attenuation and requires space- 
based observations. A number of recent detections of hot (Tgas > 
200 K) water by the Spitzer Space Telescope has demon- 
strated that abundant water vapor is a ubiquitous component in 



protoplanetary disks (jSalyk et al.||2008[ |Carr & Najita||2()08 
|Pontoppidan et al.||2010[ l. However, this emission likely arises 
within the snow-line of these young systems and therefore does 
not provide a complete picture of the distribution of water (both 
ice and vapor) in disks with sizes in excess of 100 AU. In partic- 
ular, these observations do not probe the cold (Tdust ~ 7"gas ~ 20 
K) outer parts of the disk where most of the mass resides. 

The Herschel Space Observatory (Pilbratt et al. 2010) offers 
a new opportunity for a characterization of the distribution and 
evolution of water vapor in protoplanetary disks. We report here 



a sensitive search for cold water emission in the ground-state 
lines of ortho- (o) and para- (p) H2O using the high spectral 
resolution of the HIFI instrument (d e Graauw et al.||201 0') to- 
wards the well studied DM Tau protoplanetary disk. This study is 
pait of the guaranteed time key program 'Water in Star Forming 
Regions' (van Dishoeck et al. in prep.). Stringent limits to the 
strength of both lines of a few mK suggest that the outer regions 
of the disk contain little water vapor or water ice. In Sect. 2 we 
outline the observations. Section 3 presents results from detailed 
modeling; Sect. 4 summarizes the implications of our result. 

2. Observations and results 

DM Tau is a T Tauri star located at a(2000) = 4''33"'48!7 and 
^(2000) = 18° 10' 10" with a disk diameter, estimated from CO 
emission, of ~ 1800 AU at a systemic velocity of 6 km s ' and 
i ~ -35° (jPietu et al. 2007| l. The source is a single Ml star 
dWhite & Ghez||2001| l wit h L=0.25 Lq DM T au has a chemi- 
cally rich molecular disk ( Dutrey et al.||1997 l. Accretion from 
the disk to the star also provides a source of excess UV luminos- 
ity, with an overall UV field strength of Gq ~ 240 (relative to 
the standard interstellar radiation field, ISRF;|Bergin et al. 2003 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with partic- 
ipation important from NASA. 



2004; Habing 1968). The object is a transition disk with an in- 
ner hole on the order of a few AU, based on models of Spitzer 
spectra ( |Calvet et al.|2005| l. 

DM Tau was observed with the HIFI instrument using the 
double beam switch observing mode with a throw of 3fO. On 
2010 March 22 spectra were taken in receiver band lb with an 
on-source integration time of 198 minutes, and rsys=78-96 K. 
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On 2010 March 4 spectra were taken in receiver band 4b and an 
on-source integration time of 328 minutes, and rsys=370^10 K. 

The HIFI beam of 39" at 556 GHz and 21" at 1 1 13 GHz is 
larger than the DM Tau disk with a diameter of 1277 at 140 pc. 
The beams are also larger than the pointing accuracy of Herschel 
of ~ 2". The data were recorded with Wide-Band Spectrometer 
(WBS) covering 4.4 GHz with 1.1 MHz resolution (0.59 and 
0.30 km s ' at 556 and 1113 GHz, respectively), and the High- 
Resolution Specti-ometer (HRS) covering 230 MHz at 0.25 MHz 
resolution (0. 13 and 0.067 km s ' at 556 and 1113 GHz, respec- 
tively). Both H- and V-polarizations were measured. 

The raw data were calibrated onto the T* scale by the in-orbit 
system and converted to T^\, assuming a beam efficiency of 0.74 
The data were reduced using HIPE v3.0. Subsequently, the data 
were exported to CLASSfH The HIFI flux calibration is accurate 
to 10%, while the velocity scale of HIPE v3.0 is accurate to sev- 
eral m s"'. For both lines, the WBS data were rebinned to 0.54 
km s"^ channels (or 1.0 MHz, close to the instrumental resolu- 
tion of 1.1 MHz); the HRS data were rebinned to 0.45 km s ' 
channels. All spectra, including the the H- and V-polarizations, 
were averaged together weighted by their respective noise levels. 
The resulting rms noise levels are 2.9 mK (HRS) and 1.4 mK 
(WBS) for the HjO loi-loi line, and 7.2 mK (HRS) and 4.3 mK 
(WBS) for the H2O ln-Ooo line|^ Figure 1 illustrates for band 
lb that the noise in our data decreases as (time)""-'' up to the full 
achieved integration times. 

Figure 2 presents the HRS and WBS spectra of the two water 
transitions. No strong lines are detected; for cor nparison. Fig. 2 
also shows the '^CO 1-0 spectrum of DM Tau (Kessler-Silacci 



|2004 ; Panic et al. in prep.) showing a clear emission line with 
a width of ~2 km s ' centered on the source velocity of -1-6. 1 
km s '. In the WBS spectrum of the H2O lio-loi line a weak 
feature is present between Vlsr +0.5 and -I-IO km s"', peaking 
around -t-6.6 km s ' ; a similar feature is seen in the noisier HRS 
spectrum. With T^\, - 4.3 mK, the brightest channel lies at 3cr. 
Integrated between +0.5 and +10 km s"', the feature contains 
17.2 + 3.2 mK km s"', a Scr result. A Gaussian fit to the feature 
yields best fit parameters of Vlsr - +6.8 ± 0.7 km s"', a FWHM 
width of 5.6 ± 1.2 km s"\ an intensity r^b = 2.7 ±1.4 mK, and 
an integrated intensity of 16.0 ± 2.7 mK km s"' (6cr). 

Arguments in favor of interpreting this feature as a positive 
detection of the H2O lio-loi line include the facts that the inte- 
grated intensity constitutes a 5-6cr detection and that the feature 
peaks near the systemic velocity of 6. 1 km s"' . Against the inter- 
pretation as a positive detection is a line peak that is only 2-3cr 
and a linewidth which is twice that of the '^CO line. The width 
would suggest that the emission arises from within 10 AU. With 
respect to the latter, it is interesting that HCO^ 1-0 Un e has a 
blue wing extending over ~3 km s ' ( Dutrey et al.|l997 i. 

Neither set of arguments is clearly stronger, and we interpret 
the feature in the H2O lio-loi spectrum as a very tentative de- 
tection of water vapor in the disk of DM Tau. In the remainder 
of this Letter, we will work with an intensity of < 2.7 mK for the 
H2O lio-loi line and with an upper limit to the intensity of the 
H2O lii-Ooo line of < 12.6 mK. 

Smoothing the WBS spectra to ~ 27 km s"' resolution re- 
sults in positive detections of the continuum of DM Tau at 
556.9 GHz of 3.0 + 0.5 mK (0.6 ± 0.1 Jy) and at 1113.3 GHz 



' fhttp://www.iram.£r/IRAMFR/GILDAs| 

^ Although the channel spacing in both bands in similar (0.54 vs 0.45 
km s"'), the noise in the WBS is 1.7-2.0 times lower than the noise in 
the HRS because of the larger noise bandwidth of the WBS and a V2 
loss factor in the HRS autocorrelator. 
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Fig. 1. Noise level (rms) vs integration time of the WBS and 
HRS data of the H2O lio-loi line. The noise decreases as V^- 



of 5.7 + 0.9 mK (1.3 ± 0.2 Jy). These values are consistent with 
other continuum measurements (Dutrey et al.|1996"i. 



3. Model Predictions 

3.1. Chemistry of Water Vapor in the Cold Outer Disk 

Because the midplane temperatures are well below the evapora- 
tion temperature at densities representative of the midplane of 
>150 K ( [Eraser et al.|2001||D'Alessio et al.|[2005] l, we can ex- 
pect that beyond the snow-line water is mostly frozen on the sur- 
faces of dust grains. Therefore, molecular emission arises pre- 
dominantly from the warm disk surface that is heated by stel- 
lar irradiation ( jAikawa et al.||2002| . Beyond 1 AU the dust 
in this superheated layer is heated to <100 K (Nomura et al 



2007) and again it is anticipated that water will remain as ice. 
However, the upper layers of the disk are exposed to energetic X- 
ray and FUV radiation which provide a source for non-thermal 
desorption. Models of X-ray induced desorption suggest that X- 
rays cannot release significant H2O into the gas (Najita ~et al.| 
2001 ). A more profitable method to desorb water ice in the cold 
outer disk is via photodesorpt ion, which has a m easured yield 



(molecules/photon) of ~ 10 (Obergetal 



2009 I. This has been 



suggested as providing a basal column of water vapor in the disk 
by pominik et al.| ( |20U5] l. 

Figure |3J presents results from a detailed calculation of the 
UV radiation transfer and chemistry of a standard T Tauri disk 
with Mdisk - 0.03 Mq, j^out = 400 AU and a standard gas-to- 
dust ratio (Fo gel et al.||2010 ). A key aspect in the calculation 
of the importance of photodesorption is the formation of water 
ice on the grain surface via oxygen hydrogenation and also the 
2D propagation of UV photons. This includes Ly a radiation, 
which is important for H2O ( van Dishoeck et al. 2006 1. As can be 
seen a thin layer exists where water vapor is present in moderate 
abundance (~ 10"^ - 10"'') in the cold outer regions of the disk. 



As noted by Dominik et al. ( 2005 1 the column density of 
water vapor produced by photodesorption is independent of the 
photon flux. This can be seen by balancing formation by pho- 
todesorption with destruction by UV photodissociation, giv- 
ing a maximum water abundance (Hollenbach et alT]|2009| l: 
n(U20)Jcr,F,dA = Fuvi'^s.teA^p/iceA^m'«(H20)ice. Here Fyy = 
J F^dA is the integral of the photon flux given in Fig. 3 and 
c^A ~ o"Lyff ~ 10"'^ cm"^ (van Dishoeck et aL]|2006|l. osne is 
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Fig. 2. (a)-(d) H2O spectra obtained with HIFI. A tentative fea- 
ture is present in panel (b) for H2O 1 lo-loi detected with WBS 
between Vlsr of +0.5 and +10 km s"'. A Gaussian fit to the 
feature is shown by the smooth grey line. The HRS data on this 
line in panel (a) are consistent with this result but noisier. Only 
noise is seen in panels (c) and (d) showing the ln-Ooo line. 
(e) Comparison '^CO 1-0 spectrum of DM Tau obtained with 
OVRO ( |Kessler-Silacci|2004[ Panic et al. in prep.). 



the cross section of a given site on the grain (- na^/Nsite, with 
grain radius a=0.1 jum and A^site=10^). Np=2 is a correction for 
the fact the UV photon s only penetrate the first few monolay- 
ers ( |Oberg et aL]|2009| l, with a yield of Y ~ 2 x 10 ^ /ice is 
the fraction of water ice over the total amount of ice and A^m is 
the number of monolayers. Based on this approximation we find 
;t:(H20),„„v ~ 10"^; for scaling relations see HoUenbach et al. 



(2009 1. This could be lowered if the grains have a reduced frac- 
tion of water ice or perhaps less than a monolayer of coverage. 
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Fig. 3. Abundance of water vapor (Top) and Ly a photon flux 
(Bottom) shown as a function of both radial distance and vertical 
height. The water vapor abundance is relative to total hydrogen 
and is only plotted for points where the vertical height to radial 
distance ratio is less than 1 . The line in both plots refers to the 
T = 1 surface for stellar radiation and e - I refers to the fact 
that these models were run with a standard dust to gas ratio. All 
model results taken from |Fogel et al.| ( [2010[ ). 

3.2. Comparison to Observations 

Using the above chemical model calculation as input, we use the 
molecular excitation and radiative transfer code LIME (Brinch et 
al. in prep.) to calculate the line intensity in both observed water 
lines. We use the collision rates of water with P-H2 from Faure 
et al. (2007) as provided by the LAMDA database ( [Schoier et al. 
2005 PI and convolved the results with the appropriate Herschel 
beams. In our models we assume an intrinsic broadening of 0.4 
km/s on top of a Keplerian velocity profile. 

Because of the appreciable abundance of H2O in the model 
(disk-averaged column density of 6 x lO''* cm"^), it is not sur- 
prising that significant line intensities are predicted of T^b - 
140 mK for the lio-loi line and 300 mK for the ln-Ooo line 
with a Gaussian spectral profile. Clearly, our observations rale 
out the presence of the amounts of water vapor predicted by pho- 
todesorption regardless of the details of our model. 

Absorption by low-excitation water from foreground mate- 
rial cannot explain the absence of de tected emission: the cloud 
seen in '^CO by (Dutr ey et al.|l997 i has narrow emission cen- 
tered 3 km s ' away from the DM Tau disk. Only when the col- 
umn density of water is scaled down by a factor of 130 to a 
disk-averaged value of 5 x 10'^ cm"^ does the predicted strength 
of the lio-loi line becomes consistent with the observed limits; 
the limits on the 1 1 i-Ooo line are less strict because of the higher 
noise of these observations. 

Our model could predict lines that are too intense if we over- 
estimate the collisional excitation of water. |Dick et aT] ( |2010| l 
suggested that existing collisional excitation rates for water are 
overestimated at temperatures below ~50-80 K. Decreasing the 
collisional excitation rates has little effect on the line strengths 
for the original column density. In that case, the lines are still 
highly optically thick, with maximum optical depth of 2000, and 
line trapping effectively excites the line. However, for collision 
rates lower by a factor of 10 compared to the adopted rates, re- 
ducing the disk averaged column density by a factor of 20 to 
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is sufficient to comply with the observational con- 



3x10' 
straints. 

Our generic model disk contains 0.03 Mq similar to DM Tau, 
but its 400 AU radius is only half that of DM Tau. The increased 
beam dilution only strengthens our conclusions. For a DM Tau 
specific model, Dominik et al. ( 2005 [ l predict the lio-loi line 
to be in absorption. Line profiles with combinations of emission 
and strong absorption naturally arise in disk models with strong 
temperature, density and water abundance gradients viewed un- 



der non-zero inclinations. An example is provided by Cernicharo 
[eraL] ( |2009l l for the HD 97048 disk. These models have water 
abundances peaking at much larger height above the midplane, 
explaining the low excitation resulting in strong absorption. For 
our discussion here, we stress that all these models predict emis- 
sion or absorption lines that are inconsistent with our observa- 
tional limits. Future work will focus on more detailed models 
for DM Tau, and also explore the effect of dust settling and non- 
standard gas-to-dust ratios (Brinch et al. in prep.). 

Finally, we note that Ceccarelli et al. ( 2005| l have claimed a 
4cr detection HDO in the DM Tau. This detection has been cast 
into doubt based on line formation considerations ( |Guilloteau| 
et al.|200"6l). Our models predict water will be in emission, and 



not absorption. At face value this would argue against the reality 
of the absorption. If we derive a column density from the HDO 
observations and the limits here, the D/H ratio would also be 
exceedingly high > 0.2. Deeper HDO observations are needed 
to settle this issue. 

4. Implications 

We have presented the results from a deep search for the ground 
state emission lines of 0-H2O and P-H2O towards the DM Tau 
disk. Based on the best theoretical knowledge we have to date, 
water vapor should be present in the outer disk and presum- 
ably emissive. However, our sensitive observations show that, 
at least for this object, it is not. Our limit on the P-H2O line pre- 
cludes an extreme o/p ratio as explanation for the low 0-H2O 
emission strength. There are two potential explanations for this 
result. Either water vapor is unemissive as would be the case 
if the excitation at low temperature is lower than generally as- 
sumed ( IDick et al.|2010| ). [Cernicharo et al.| ( |2009 ) also show that 
the water excitation sensitively depends on the adopted collision 
rates with o- and /7-H2. However, our calculations suggest that 
line trapping is sufficiently effective at the predicted water abun- 
dances to still produce detectable lines. Alternatively, our physi- 
cal/chemical understanding may be incorrect. The photodesorp- 
tion yield is measured in the laboratory at low temperature and 
the abundance is independent of the photon flux so these aspects 
appear unlikely to provide the answer. One intriguing possibility 
is that the upper layers of the outer disk are 'dry' which could 
be the case if only bare grains are present in the region where 
UV photons are present. A well known key aspect of disk phys- 
ical evolution is the coagulation and settling of dust grains to 
the disk midplane ([Weidenschilfing & Cuzzi 1993[ DuUemond 



|& Dominik|2004[ [Furlan et al.|2006| ). Icy grains present a more 
favorable surface for grain coagulation and would therefore be- 
come larger and settle to deeper layers than their bare silicate 
counterparts (Dominik & Tielens|1997| l. While upward mixing 
of gas and small grains may occur, larger ice-bearing grains re- 
main in the midplane. In addition, the total ngnin CTgr^in will be 
reduced, thereby reducing the efficiency of grain surface forma- 
tion of H2O upon which the photodesorption model depends for 
water vapor creation (Hollenbach et al. 2009 ). This 'cold-finger' 
effect was also proposed by |Meijerink et al.| ( |2009j ) to explain the 



truncation of warm water vapor beyond ~ 1 AU seen in Spitzer 
measurements. Thus, Herschel and Spitzer both suggest that the 
disk around DM Tau is settled. In summary, our Herschel results 
suggest that less than 1-5% of the water ice reservoir survives 
in the UV-illuminated outer disk regions around DM Tau. If this 
finding is confirmed by more detailed models and by additional 
observations, Herschel may be telling us something entirely new 
about the chemical structure of protoplanetary disks. 
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